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Abstract 

The well known Sunyaev-Zel'dovich (SZ) effect is reexamined using a 
Doppler shift type mechanism arising from the scattering of photons by 
electrons in an optically thin gas. The results are in excellent agreement 
with the observational data as well as with the results obtained using the 
diffusive pictures. A comparison of the results here obtained with other 
approaches is thoroughly discussed, as well as some important extensions 
of this method to other aspects of the SZ effect. 



1 Introduction 

As it is well known today, the cosmic microwave background radiation (CMBR) 
that fills the universe exhibits an almost perfect blackbody Planckian spectrum 
corresponding to a temperature of 2.726 K. A little over thirty years ago, it was 
predicted by the astrophysicists R.A Sunyaev and Ya. B. Zel'dovich that this 
spectrum could be distorted when photons from this radiation would penetrate 
large structures such as the hot intracluster gas now known to exist in the uni- 
verse [1). This distortion would arise from the interaction between the photons 
and electrons which constitute the plasma responsible for an important main 
contribution to the total mass contained in the cluster. Such distortion is only 
a very small effect changing the brightness of the spectrum by a figure of the 
order of 0.1 percent. This effect is now called the Sunyaev-Zel'dovich effect and 
its detection is at present a relatively feasible task due to the modern observa- 
tional techniques available. Its main interest lies on the fact that it provides 
information to determine important cosmological parameters such as Hubble's 



1 



constant and the baryonic density. Since all these and other important facts are 
readily available in the somewhat broad literature on the subject Q-Q, we will 
not pursue any further details of its cosmological implications. 

In this work we are mainly concerned with the nature of the various inter- 
pretations that have been provided to explain the effect and to propose another 
one which we believe, is much simpler to grasp than the others, and offers a 
much more direct way to account for the relativistic corrections, as well as the 
influence of z, the redshift factor, in the equations for the distorted spectrum. 

The main idea behind the method we want to discuss is that photons incident 
on the hot electron gas of the plasma change their frequency by an absortion- 
emission process, so that the intensity of the spectral line is changed through 
the Doppler effect. 

The overall effect is due to those photons that happen to be captured by an 
electron, an event which is in general unlikely to occur specially in an optically 
thin gas. Thus, an electron moving with a given thermal velocity emits (scatters) 
a photon with a certain incoming frequency Vo and outgoing frequency u. The 
line breath of this process is readily calculated from elementary kinetic theory 
taking into account that the media in which the process takes place has an 
optical depth directly related to the usual Compton parameter y. When the 
resulting expression is convoluted with the incoming flux of photons obtained 
from Planck's distribution, one easily obtains the disturbed spectra. The results 
derived from the final equation are shown to be in agreement with those obtained 
from the observational data. To present these results we have divided the paper 
as follows. In section II a brief review of the previous methods developed to 
explain the SZ effect is given. This will allow a thorough comparison with our 
proposal. In Section III we develop our ideas and give the comparison between 
the mathematical results with both, those derived before and the observational 
ones. Section IV is left for some concluding remarks and future directions of 
this work. 



2 The Diffusion- Scattering Picture of the SZ ef- 
fect. 

As it was clearly emphasized by the authors of this discovery in their early 
publications ||^ as well as by other authors, the distortion in the CMBR 
spectrum by the interaction of the photons with the electrons in the hot plasma 
filling the intergalactic space is due to the diffusion of the photons in the plasma 
which, when colliding with the isotropic distribution of a non-relativistic electron 
gas, generates a random walk. The kinetic equation used to describe this process 
was one first derived by Kompaneets back in 1956 0-[||. For the particular 
case of interest here, when the electron temperature Te ^ IQ^K is much larger 
than that of the radiation, T = Tj^ad (2.726 K), the kinetic equation reads |^ 

dN _^^cPN ^^dN 
dy dv'^ dv 
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Here N is the Bose factor, = (e^ — 1) ^ , x — j^^, v is the frequency, h 
is Planck's constant, k Boltzmann's constant and y the "Compton parameter" 
given by, 

kTe f kTe , , 

y = 2 / "eC at = -T (2) 



me being the electron mass, c the velocity of light, ctt the Thomson's scattering 
cross section and Ue the electron number density. The integral in Eq. is 
usually referred to as the "optical depth, r" , measuring essentially how far in 
the plasma can a photon travel before being captured (scattered) by an electron. 
Without stressing the important consequences of Eq. readily available in 
many books and articles we only want to state here, for the sake of 

future comparison, that since the observed value of No{v) is almost the same as 
its equilibrium value Neq{v), to a first approximation one can easily show that 



(e^ - 1) 
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(3) 



implying that, in the Ray leigh- Jeans region {x << 1), 



SN 

and in the Wicn limit {x >> 1) : 

If we now call loii^) the corresponding radiation flux for frequency v, defined as 

^o(^) = -^U.iT) (6) 

where Uy{T) = ^"^^Z Nf.q{v) is the energy density for frequency v and temper- 
ature T , and noticing that 



5T _ (d{\nl{v))\ (51 



T V 9(lnT) J \ I 



(7) 



we have that the change in the background brightness temperature is given, in 
the two limits, by 

— ^-2y,x«l (8) 

and 

ST 

— =xy,x»l (9) 

showing a decrease in the low frequency limit and an increase in the high fre- 
quency one. Finally, we remind the reader that the curves extracted from Eq. 
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(H) for reasonable values of the parameter "y" show a very good agreement with 
the observational data. 

Nevertheless, many authors, including Sunyaev and Zel'dovich themselves, 
were very reluctant in accepting a diffusive mechanism as the underlying phe- 
nomena responsible for the spectrum distortion. The difhculties of using diffu- 
sion mechanisms to study the migration of photons in turbid media, specially 
thin media, have been thoroughly underlined in the literature Sev- 
eral alternatives were discussed in a review article in 1980 |^ and a model was 
set forth by Sunyaev in the same year |l^ based on the idea that Compton 
scattering between photons and electrons induce a change in their frequency 
through the Doppler effect. Why this line of thought has not been pursued, 
or at least, not widely recognized, is hard to understand. Two years ago, one 
of us (ASV) [|4| reconsidered the single scattering approach to study the SZ 
effect. The central idea in that paper is that in a dilute gas, the scattering law 
is given by what in statistical physics is known as the dynamic structure factor, 
denoted by S{k, v) where k = and A is the wavelength. In such a system, 

this turns out to be proportional to exp(— where w is the broadening of the 
spectral line given by 

^^2^^)1/2^ (10) 

C TO 

If one computes the distorted spectrum through the convolution integral 

/■oo 

I{J^)= Io{D)S{k,i>~{l-ay)iy)di? (11) 







where the corresponding frequency shift ^ has been introduced through Eqs. 
(H^) (a — —2 + x), one gets a good agreement with the observational data. This 
is exhibited in Fig. (1) for y = 10~^. This result is interesting from, at least, 
two facts. One, that such a simple procedure is in agreement with the diffusive 
picture. This poses interesting mathematical questions which will be analyzed 
elsewhere, specially since the exact solution to Eq. (1) is known (see Eq. (A-8) 
ref. 1^). The other one arises from the fact that this is what triggered the idea 
of reanalyzing the SZ effect using elementary arguments of statistical mechanics 
and constitutes the core of this paper to be presented next. 



3 The Doppler Effect Approach. 

We begin this section by simply reminding the reader that if an atom in an ideal 
gas moving say with speed in the x direction emits light of frequency Vo at 
some initial speed Ux{0), the intensity of the spectral line / (v) is given by 



Io{u) 



exp 



2kT 



(12) 
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Figure 1: Comparison of the CMBR SZ distortion as computed by Eqs. (3,6) 
(solid line) and Eq. (11) (dashed line) i.e using the approach of Ref. [13] , with 
y = 10^'"' and T — lif'K. SI is measured in erg cm~^ ster~^, notice that 
the figure is scaled by a factor of 10^*. 



Eq. (g_2|) follows directly from the fact the velocity distribution function 
an ideal gas is Maxwellian and that Ux, Ux{0) and are related through the 
Doppler effect. For the case of a beam of photons of intensity lo {i') incident 
on a hot electron gas regarded as an ideal gas in equilibrium at a temperature 
Te the full distorted spectrum may be computed from the convolution integral 
given by 



1 



lo [v) exp 



/(y)^ 



W{v) 



dv 



(13) 



which defines the joint probability of finding an electron scattering a photon 
with incoming frequency D and outgoing frequency v multiplied by the total 
number of incoming photons with frequency v. (v) is the normalizing 

factor of the Gaussian function for f{y) — 1, W (v) is the width of the spectral 
line at frequency v and its squared value follows from Eq.(hd) 



(i^) 



4:kT, 



(14) 



where r is the optical depth whose presence in Eq. ( p^ ) will be discussed later. 
The function /(y) multiplying v in Eq. ( p^ is given by f{y) = 1 + ay where 
a = —2 in the Rayleigh-Jeans limit and a = xy in the Wien's limit, according to 
Eqs. (8-9) and the fact that ^ = ^ for photons. Eq.(|l^) is the central object 
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of this paper so it deserves a rather detailed examination. In the first place it 
is worth noticing that Io{v), the incoming flux, is defined in Eq. (|6|). Secondly, 
it is important to examine the behavior of the full distorted spectrum in both 
the short and high frequency limits. In the low frequency limit, the Rayleigh- 
Jeans limit Io{y) = '^^'^i , so that performing the integration with a — —2 and 
noticing that y is a very small number, one arrives at the result 

( p^ is in complete agreement with the value obtained using Eq. ( ^ , the photon 
diffusion equation. In the high frequency limit where a = xy and Io{v) = 
^^5^e~^, a slightly more tedious sequence of integrations leads also to a result 
at grips with the diffusion equation, namely 

— = x'^y, x»l (16) 

Why both asymptotic results, the ones obtained with the diffusion equation and 
those obtained from Eq. (^3|) agree so well, still puzzles us. At this moment 
we will simply think of them as a mathematical coincidence. Nevertheless, it 
should be stressed that in his 1980 paper, Sunyaev reached rather similar 
conclusions although with a much more sophisticated method, and less numeri- 
cal accuracy for the full distortion curves. The distorted spectrum for the CMBR 
radiation may be easily obtained by numerical integration of Eq. ( p^ ) once the 
optical depth is fixed, y is determined through Eq.(||) and a = —2 + x. The 
intergalactic gas cloud in clusters of galaxies has an optical depth r ~ 10^^ ||l^ . 
To compare our results with the curves obtained from the traditional approach, 
we have plotted 61 (f) for several temperatures, typical of the hot intra-cluster 
gas, in the non-relativistic range, in Figs. (2) and (3). 



From the results obtained one appreciates the rather encouraging agreement 
between the observational data and the theoretical results obtained with the 
three methods, the diffusion equation, the structure factor or scattering law 
approach, and the Doppler effect. This, in our opinion is rather rewarding and 
some efforts are in progress to prove the mathematical equivalence of the three 
approaches. From the physical point of view, and for reasons already given 
by many authors, we believe that the scattering Doppler effect picture does 
correspond more with reality, specially for reasons that will become clear in the 
last section. 



4 Discussion of the results 

Eq. (p^, the main result of this work hardly needs a more detailed explanation, 
it is a direct result of elementary statistical mechanics, except of course for 
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Figure 2: Comparison of the CMB SZ distortion as computed by Eqs. (3-6) 
(dashed hue) and Eq. (13) (solid hue), with r = lO'^ , T = 2KeV (lower 
curves), T = 3KeV (middle curves) and T = AKeV (upper curves). Here 
a — —2 in Eq. (13). Since these are actually six curves, it is clear that the 
accuracy achicvcid by Eq. (13) is remarkable. SI is measured in erg cm~'^ 
ster~^, the figiuc is again scaled by a factor of 10"^^. 
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Figure 3: The same as in figure 2, but taking r = 7 x 10 ^ . The respective 
temperatures are 5, and 7 KeV. 



7 



the factor r which appears in the broadening of the spectral hne. That this 
broadening must somehow depend on r is clear. As argued by many authors 
in the case of clusters of galaxies, the plasma is optically thin, so most of the 
photons are not scattered, and how many of them are must depend on the 
optical depth. Now, why in particular Eq. (|l^) is valid remains to be rigorously 
shown. In our case it came as a mere accident since in ref. the mass taken 
in w to perform the calculations was that of a proton, which turns out to be 
of the same order of magnitude as ^ for y = 10~^. There is in fact one way 
of understanding this puzzle. In a pure diffusive (Fickian) process described 
by a Gaussian function one knows that the squared value of a single line width 
grows as t. Looking at the original equation of Kompaneets and thinking only in 
its diffusive terms, one sees that the variable time is replaced by the Compton 
parameter "y", so that one should expect that, to a first approximation, the 
width of the curve grows as ^/y. Hence, inside the gas, the effective frequency 
at which the photon propagates through the gas is, according to Eq. (p^, 
proportional to x v. This argument, which is based on ideas pertinent to 
the diffusion mechanism must be extended to the Doppler picture, but this has 
not yet been accomplished. 

The two main advantages of Eq. ( p^ ) are, firstly, that its generalization 
to the relativistic case is straightforward and indeed, since in relativity theory 
one has both the transverse and parallel Doppler effects one can easily 
study what their influence is, if any, in the full distorted spectrum. For the 
parallel effect the calculations carried out so far agrees reasonably well with 
those reported by Rephaeli and other authors Q and will be the subject of a 
forthcoming publication. The case of the transverse effect is under study. Sec- 
ondly, in all the work we have reported here, the cluster of galaxies, the source 
of the scattering processes of the CMBR and the electrons, has been consid- 
ered a system at rest. One can now extend the calculations assuming that the 
cluster recedes to infinity according to Bubble's law and introduce the redshift 
factor into the formalism. The resulting prediction could be compared with 
the observational data and extract valuable information on such an important 
cosmological parameter. Work in this direction is also in progress. Concluding, 
we have the strong conviction that the Doppler effect approach to study the 
thermal SZ effect is not only very clear from the physical point of view, it also 
contains ingredients which are promising for further studies more advantageous 
than those offered by the diffusive model. 

References 

[1] Y.B. Zel'dovich and R.A. Sunyaev, Astrophys. Space Sci. 4 , 301 (1969). 

[2] R.A. Sunyaev and Ya. B. Zel'dovich. Comm. Astrophys. Space Phys. 4, 173 
(1972). 

[3] P.J.E. Peebles, Physical Cosmology (Princeton Univ. Press, Princeton, 2nd. 
Ed. 1993). 



8 



[4] MM. Joy and J.E. Carlstrom; Science 291, 1715 (2001) and ref. therein. 

[5] R.A. Sunyaev and Ya. B. Zel'dovich; Ann. Rev. Astroni. Astrophys. 18, 537 
(1980) 

[6] S.H. Hansen, S. P astor and D.V. Semikoz, Astrophys. J. 573, (2002) L69 
iastro-ph/0205295| 

[7] A.S. Kompaneets; JETP. 4, 730 (1957). 

[8] R. Weymann; Phys. Fluids 8, 2112 (1965) 

[9] D.J. Durian and J. Rudnick; J. Opt. Soc. Am. , A14 , 235 (1997) 

[10] A. H. Gandibakhche and G. H. Weiss in Progress in Optics, XXXIV; E. 
Wolf, Ed. (Elseveier, Amsterdam 1995), pp 332. 

[11] G. H. Weiss; an Overview of Persistent Random Walks and the Telegra- 
pher's Equation in Physycs, NIH- Washington, preprint (2001). 

[12] R.A. Sunyaev; Sov. Astron. Lett. 6, 212 (1980) 

[13] Y. Rephaeh, Astrophys. J. 445, 33 (1995); 

A.D. Challinor and A.N. Lasenby, Astrophys. J. 499, 1 (1998); 
N. Itoh, Y. Kohyama, and S. Nozawa, Astrophys. J. 502, 7 (1998); 
S. Nozawa, N. Itoh, and Y. Kohyama, Astrophys. J. 507, 530 (1998); 
A.D. Dolgov, S.H. Hansen, D.V. Semikoz and S. Pastor, Astrophys. J. 554 
74 (2001) ]astro-ph/00104T^ . 

[14] A. Sandoval- Villalbazo; Physica (2002) (in press). | |cond-mat/0109274| ]. 

[15] C. Mdller; The theory of Relativity (Oxford Univ. Press. Oxford, 1952). 



9 



